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Abstract 

Directed and elliptic flow of inclusive photons near mid-rapidity in 158 A GeV 
Pb+Pb collisions has been studied. The data have been obtained with the pho- 
ton spectrometer LEDA of the WA98 experiment at the CERN SPS. The flow 
strength has been measured for various centralities as a function of px and rapidity- 
over 0.18 < pr < 1.5GeV/c and 2.3 < y < 2.9. The angular anisotropy has been 
studied relative to an event plane obtained in the target fragmentation region that 
shows the elliptic flow to be in-plane. The elliptic flow has also been studied using 
two-particle correlations and shown to give similar results. A small directed flow 
component is observed. Both the directed and elliptic flow strengths increase with 
Pt- The photon flow results are used to estimate the corresponding neutral pion 
flow. 

Key words: ultrarelativistic heavy-ion collisions, collective flow, inclusive photons 



3 



1 Introduction 



Heavy ion collisions at relativistic energies provide a means to study the prop- 
erties of nuclear matter at high temperature and density. In such collisions 
it is expected that a high density interaction zone is formed. If this system 
thermalizes, the thermal pressure will generate collective transverse expansion 
[1]. Such collective flow, and especially its anisotropy, will reflect the time evo- 
lution of the pressure gradients of the system and can provide information on 
the equation of state (EOS) in the initial phase [2,3] and during the expansion 
[4], and in particular about the possible formation of the Quark Gluon Plasma 
(QGP) [5,6]. 

QGP formation requires high energy density and local thermalisation of the 
system. While transverse energy measurements indicate the attainment of high 
energy densities [7] , the degree of thermalization has not been unambiguously 
determined. Measurements of collective flow may be one of the strongest hints 
related to the degree of thermalisation. One may consider two extreme cases: 
the low density limit, where the mean free path is comparable or larger than 
the system size for which cascade models may be appropriate, and the hydro- 
dynamic limit, where the mean free path of the particles is much less than 
the system size. The study of collective flow as a function of beam energy and 
system size, transverse and longitudinal momentum, and particle species may 
allow to separate these two scenarios and to recognize the hadronic or partonic 
(QGP) nature of the reaction [8]. 

The anisotropic flow of charged fragments has been measured in nuclear colli- 
sions at various beam energies: at 0.1-1.0 A GeV [9,10,11], at 10A GeV [12], at 
158A GeV [13,14,15,16,17] and also at ^IJ^ = 130 and 200 GeV [18,19,20]. 
The first j>r-integrated photon flow measurement was performed by the WA93 
experiment [15] for 200 A GeV S+Au collisions, where photons were measured 
with the Photon Multiplicity Detector (PMD). In the present work we report 
first results of photon collective flow measurements in 158 A GeV Pb+Pb 
collisions using photons identified and momentum analyzed in the lead glass 
calorimeter LED A of the WA98 experiment. Preliminary results have been 
presented in [21,22]. 



2 Experimental Setup 

The data presented here were obtained in the WA98 experiment [23] for 158 
A GeV Pb+Pb collisions at the CERN SPS. The WA98 setup consisted of 
large acceptance hadron and photon spectrometers, calorimeters for forward 
and transverse energy measurements, and detectors for photon and charged 
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particle multiplicity measurements. 



The centrality of the event was determined by the total transverse energy, 
Et, measured with the mid-rapidity calorimeter (MIRAC) [24], which covered 
the pseudo-rapidity range of 3.5 < 77 < 5.5. It was placed at 24.7 meters 
downstream from the target and consisted of 30 stacks, each divided vertically 
into 6 towers, of size 20 x 20 cm 2 each. The MIRAC measured both the 
transverse electromagnetic E T m and hadronic E T ad energies. Events with large 
Et correspond to the most central collisions with small impact parameter. 
The minimum bias trigger required a beam trigger with a MIRAC transverse 
energy greater than a low threshold. 

The Plastic Ball spectrometer, used for the reaction plane (RP) determination, 
had full azimuthal coverage in the pseudorapidity range of —1.7 < r\ < 0.5 
(i.e. in the target fragmentation region with polar angles 70° < 9 < 160°). 
It consisted of 655 detector modules and allowed to identify pions, protons, 
deutrons, tritons, 3 He, and 4 He with kinetic energies of 50 to 250 MeV by 
the AE — E method. Each module comprised a slow 4 mm thick CaF 2 AE 
scintillator followed by a fast plastic scintillator, both read out by a common 
photomultiplier [25]. 

Photons, of which ~ 85% originate from 7r° decay [26], were detected in the 
electromagnetic calorimeter LED A, a highly segmented photon detector lo- 
cated 22.1 m downstream of the target and covering the photon rapidity re- 
gion 2.3 < y < 2.9, i.e. backwards of mid-rapidity (y = 2.9). LEDA consisted 
of 10080 TF1 lead-glass modules, read out by FEU-84 photomultipliers. The 
photomultiplier high voltage was generated on-base with custom developed 
[27] Cockcroft- Walton voltage-multiplier type bases which were individually 
controlled by a VME processor. The photomultiplier signals were digitized 
with a custom-built ADC system [28]. The dimensions of each module were 
4x4x40 cm 3 (14.3 radiation lengths depth and 1.1 Moliere radius width). Each 
group of 24 modules had its own calibration and gain monitoring system based 
on a set of 3 LEDs mounted inside a sealed reflecting cover dome. Each mod- 
ule viewed the reflected LED light through an aperture on the front surface, 
while the LED light was simultaneously monitored by a PIN-photodiode [29]. 

LEDA was calibrated with 10 GeV electrons in the XI beam at the CERN SPS 
in the years 1993-1994. Electron beams with energies from 3 GeV to 20 GeV 
were used to measure the energy and position resolution, and the energy non- 
linearity. The measured energy resolution was a/E = (5.5±0.6)%/a/E + (0.8± 
0.2)% and the measured position resolution was a/E = (8.35±0.25)mm/v^E'+ 
(0.15 ± 0.07)mm. A more detailed description of the WA98 setup is given in 
[16,23,26]. 
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3 The Methods 



Two complementary methods have been used in this analysis of anisotropic 
flow: the study of single-particle angular distributions with respect to an es- 
timated reaction plane (reaction plane method) and the study of two-particle 
correlations (correlation method). Both methods assume that the underlying 
azimuthal distribution of particles with respect to the reaction plane (RP), 
which is the plane that contains the impact parameter and beam direction 
vectors, can be described by a Fourier decomposition: 



where is the azimuthal angle of the emitted particle, and is the az- 
imuthal angle of the RP. The anisotropic flow is characterized by the values 
of the Fourier coefficients v± : for directed flow, and v 2 , for elliptic flow. 

3.1 The Reaction Plane Method 

The conventional reaction plane method [30,31] uses the distribution of parti- 
cles in their azimuthal angle relative to the estimated reaction plane. Because 
the true reaction plane is not known in the experiment, one has to establish an 
event plane (EP) from the measured particles as an estimate for the reaction 
plane. 

Particles measured in the target fragmentation region show significant directed 
flow - proton, deuterons, and heavier fragments are emitted in the reaction 
plane in one direction, while pions are emitted in the opposite direction, as 
has been measured with the Plastic Ball detector [32]. This information has 
been used to calculate an event plane angle from the particles measured in 
the Plastic Ball: 



where the sum runs over all fragments and identified positive pions. E l T and 
(pi are the transverse kinetic energy and azimuthal angle in the laboratory 
frame of the i— th particle, respectively. For pions, was replaced by + n to 
account for the opposite sign pion directed flow [32]. 

The azimuthal distribution of photons detected by the LEDA calorimeter has 
been studied relative to this event plane angle. The distributions are studied 



1 dN 



1 + 2ui cos (0 - y RP ) + 2v 2 cos (2(0 - V RP )) , 



(1) 




(2) 
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as a function of A$ = 7 — & E p and are fitted with: 
1 dN 



NdA<5> 



= 1 + 2< s cos (A$) + 2< s cos (2A$) . (3) 



The measured EP doesn't coincide exactly with the true RP because of the 
finite number of detected particles and resulting fluctuations. Because of this 
finite reaction plane resolution, the coefficients obtained from the fits have to 
be corrected by dividing them by the event plane resolution correction factors 
(RCF n ): 



yobs 



RCF, 



(4) 



n 



The resolution correction functions RCF n are given by [31]: 



RCF n = (co S (n($ EP -V RP ))) 



2^/2 



Xmexp 





(5) 



where ^ rp is the true RP angle, Xm is the resolution parameter, which is 
proportional to the square root of the multiplicity, m is the order of the Fourier 
component used for calculation of the event plane and k — n/m. 

Since ^ rp is unknown, the RCF's must be determined from the measured EP 
themselves. This can be done by a subevent analysis in which each event, in 
this case consisting of hits in the Plastic Ball detector, is randomly divided 
into two subevents (A and B) and for each subevent the EP angle $ep ($^4 or 
$b) is calculated. The quantity (cos(n($ 7 4 — $b))) is determined directly from 
the subevent correlation function. It is then used in equation 5 to obtain the 
parameter Xm b f° r the subevent multiplicity, which is then used to calculate 
Xm = V^Xm b f° r the full event, and finally to obtain RCF n . In this analysis 
we have used m — 1 (k — 1) for RCF ± and for RCF 2 (k = 2). 



3.2 The Correlation Method 



Alternatively, the flow values have been obtained from the azimuthal correla- 
tions of photon pairs as described in [2,31,33,34,35]. The correlation function 
is calculated as: 

_ d 2 N/d(j) 1 d<p 2 

C " m = dN/d^-dN/d^ (6) 
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where A0 = <f>i— 2 . The C 77 are calculated as the ratio of the true two-photon 
distribution to the pair distribution from mixed events. This is necessary to 
correct for distortions from the limited acceptance of the photon detector. In 
the event-mixing procedure we have taken care to mix only events with similar 
global properties and to use identical cuts, especially regarding the two-cluster 
separation within the detector. The two-photon correlation functions have 
then been fitted with: 

C 77 (A0) = 1 + 2v\ cos(A0) + 1v\ cos(2(A0)). (7) 

From these fits, no significant v\ component was observed. Since very little 
directed flow is expected near mid-rapidity, we have set v\ = in this analysis. 

If collective flow is dominant, this method should be equivalent to the previous 
one since (a) the correlation between every particle and the RP induces a 
correlation amongst the particles, and (b) correlating two subevents amounts 
to summing two-particle correlations [36]. The correlation method has the 
advantage that no EP determination, and therefore no resolution correction, 
is needed. However, non-flow correlations, such as back-to-back correlations 
due to momentum conservation [22], should be taken into account. 



4 Analysis 

The data presented here were taken in 1995 and 1996 with the 158 A GeV 
Pb ion beams of the CERN SPS. Pb targets of 495 and 239 mg/cm 2 thickness 
were used. About 10 7 events were analyzed. The events were divided into 8 
centrality classes defined by intervals in the total E T measured by MIRAC, as 
summarized in Table 1. The centralities are expressed as fractions of the min- 
imum bias cross sections as a function of the total Et, measured by MIRAC. 
In addition, a multiplicity of greater than 3 fragments measured in the Plastic 
Ball was demanded to allow for a reasonable determination of the EP. This 
has a significant effect for the peripheral bins and causes a slight bias towards 
higher multiplicity within the bin. 

In addition to the 8 classes of centrality shown in Table 1, studies were per- 
formed with combined centrality classes: 2 + 3 (47-83%), 4 + 5 + 6 (13-47%) 
and 7 + 8 (0-13%). The table also shows the number of participants N part as 
calculated in a Glauber type calculation discussed in [37] . 

To further suppress the hadron contamination in the photon sample, only 
those showers in the calorimeter have been used which satisfy the following 
cuts: 
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Et class 
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83 — 100% 


U — Zo.00 


i n -i-9 

1U zalL 


2 


bo — 83% 


zo.oo — iy.uo 


98 -1-9 


3 
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Oo ItZ 
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9 /I A 70/ 

Z4 — 4 / /o 


1(51 ^ 981 
101.00 — Zol.UO 


loo ZtO 


5 


19 - 24% 


981 n 1 ^ ^i 8 


9nc; _i_o 


6 


13 - 19% 


318.05 - 361.55 


247 ±2 


7 


6.5 - 13% 


361.55 - 410.95 


291 ±2 


8 


- 6.5% 


> 410.95 


351 ±1 



Table 1 

Centrality classes used in this analysis. The percentage of the measured minimum 
bias cross section included in each class is given. Also given is the corresponding 
average number of participants for each class with an estimate of the systematic 
error. Cuts on Et are for 1995 data set, 1996 data set cuts are shown in [38]. 

• The measured energy was greater than 0.75 GeV. This cut suppressed min- 
imum ionizing particles. 

• The lateral dispersion of the shower was less than a maximum value. This 
suppressed showering hadrons. 

These cuts kept the hadron contamination in the photon sample to less than 
w 7% [26,39]. 

The observed raw $£;p-distributions showed a variation due to detector biases, 
such as dead channels and inefficiency, of less than 5%. This non-uniformity 
has to be removed before extracting the flow strength from the measured 
correlation functions. This has been done by two methods: 

• The real distribution was divided by the equivalent distribution for mixed 
events, where <p 1 and fgp are taken from different events. 

• In accumulating the distributions, the entries were weighted with the inverse 
of the $ep distribution. 

Both methods gave consistent results within errors, we have used the second 
method for the final result. 

Fig. 1 (left) shows examples of the measured photon azimuthal correlation 
functions for different centralities. A clear modulation is seen, especially in 
the semi-central classes, which can be described well with fits to equation 3. 

The resolution correction factors RCF's for each centrality were determined 
as described above from the subevent correlation functions as shown in Fig. 1 
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(right). In this case the Plastic Ball acceptance correction was done by using 
mixed subevents, i.e. subevents from different events. The deviation of the 
mixed subevents from unity, i.e. the order of magnitude of this correction, was 
less than 2% (as shown in Fig. 1 (right)). 

Fig. 2a shows the values of the resolution correction factors RCF's determined 
in this way as a function of the number of participants. A stronger subevent 
correlation implies a better determination of the RP, and the corresponding 
values of the RCF's are larger [40]. It is seen, that the subevent correlation 
is strongest for semi-central events, whereas for peripheral and central events 
the quality of the EP - determination is worse. The RCF's shown in Fig. 2a 
have been obtained for the 1996 beam time. The RCF's have been calculated 
and applied independently for the 1995 beam time, with values found to be 
smaller by ~ 20 - 40%. 

The v 2 elliptic flow values for photons have also been obtained from the two- 
photon correlation functions for various centrality classes, although not all of 
these precisely match those used for the reaction plane method analysis [41]. 
The V2 elliptic flow values for photons obtained by the two methods are com- 
pared in Figs. 2b and 3. The two photon correlation functions have been fitted 
with Eq. 7 to extract i>2 with v\ = 0. Fig. 2b shows the centrality dependence 
of i> 2 integrated over p T > 0.18 GeV/c. The solid circles are the values ob- 
tained from the reaction plane method and the open circles are those from 
the correlation method. For centralities corresponding to N part > 50 there is 
very good agreement between the two methods. Fig. 3 shows the elliptic flow 
as a function of pt for central and two semi-central classes. There is good 
agreement between the methods. 

One should note that the most peripheral classes suffer from several problems: 
The determination of the reaction plane has a large uncertainty, especially 
fluctuations due to the low multiplicity play a role here. In addition, for the 
correlation method it has been seen that there may be considerable non-flow 
components for these centralities [22] which would influence the extraction 
of the flow. Furthermore, we have neglected directed flow in the correlation 
method which may also lead to systematic errors in the v-i determination. The 
results from the reaction plane method are used in the discussions that follow. 

As a consistency check we have performed the reaction plane analysis inde- 
pendently on the data from the two different beam periods, where the quality 
of the RP determination was found to be different. The resulting flow values 
are in good agreement. 
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5 Systematic Error 



The systematic errors of the obtained coefficients include: 

• uncertainties in the event plane determination in the Plastic Ball due to non- 
uniformity of acceptance and efficiency or imperfect particle identification, 
estimated as < 2%. 

• uncertainties in the photon angular distributions due to charge particle con- 
tamination of photons in LEDA. This contamination is less than ps 7% 
([26]). The uncertainties are related to a difference in the observed flow be- 
tween charged pions and pion decay photons and were estimated as < 3.5%. 

Other sources of the systematic errors have been investigated by comparing 
the results obtained under different conditions. The following checks have been 
performed: 

(1) Different weights (E T vs. pr) have been used for the determination of the 
event plane. 

(2) Different acceptance regions of the Plastic Ball detector (70° < 9 < 160° 
or 60° < 6 < 160°) have been used. 

(3) Different identification cuts for the Plastic Ball fragments were used. 

(4) The event plane was determined with and without including pions in the 
Plastic Ball. 

(5) Results of the two different beam periods were analyzed separately. 

Another possible contribution to the systematic error is due to non-flow cor- 
relations. Among such non-flow effects relevant for azimuthal correlations are 
the correlations due to momentum conservation, long- and short-range two- 
and many-particle correlations (due to quantum statistics, resonances, jet or 
mini-jet production, etc.). The contribution of non-flow correlations scales 
as 1/N, where N is the multiplicity of particles used to determine the event 
plane. The "momentum conservation" contribution increases with the fraction 
of particles detected, and the relative contribution of Bose-Einstein correla- 
tions would be independent of N. The effect of non-flow correlations in the 
reaction plane method is expected to be small as the event plane is determined 
from particle which have a large pseudorapidity separation Ai] > 1.8 from the 
photons. The non-flow contributions to the event plane determination have 
been investigated by studying the dependence of the RCF's on the Plastic 
Ball multiplicity. The conventional subevent analysis was performed exclud- 
ing various fractions of the Plastic Ball particles to investigate the deviation 
of the flow parameter Xm from Xm ~ ^/N (see eq. 5). 

All of these systematic checks lead to a total systematic error estimate of less 
than ±17% for both measured values of the flow coefficients. Unless explicitly 
stated otherwise, the systematic errors are not included in the figures. 
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6 Results 



Fig. 4a,b shows the final values of the directed flow v\ and the elliptic flow 
v 2 coefficients as a function of the number of participants. The data were 
integrated over p T > 0.18 GeV/c and over rapidity y = 2.3 — 2.9. 

Both types of flow decrease in strength with the number of participants. The 
sign of all values is positive with the sign convention that positive V\ corre- 
sponds to the directed flow direction of the protons in the projectile fragmen- 
tation direction and positive v 2 corresponds to the in-plane direction. Thus, 
the directed flow of photons below mid-rapidity is in the same direction as the 
directed flow of pions, oppositee to the protons, in the target fragmentation 
region, and the elliptic flow is oriented in the reaction plane. 

The transverse momentum dependence of the photon flow integrated over 
y = 2.3 — 2.9 is shown in Fig. 4c, d for different centrality classes. Both flow 
coefficients show an increase with p T which is compatible with a blast wave 
behavior (see below, [19], [42]). 

To be able to compare the flow of photons to other experimental results we 
have attempted to establish the relation between the observed photon flow and 
the underlying flow of the parent neutral pions. This has been investigated in 
Monte-Carlo simulations. In the simulations, 7r°'s were generated according 
to the pt spectrum measured by WA98 [26]. The 7r° azimuthal distributions 
were modulated with directed and elliptic flow components. The values of the 
7T° Vi and t>2 and their pr- dependence have been constrained by the photon 
measurement. The pr-dependence of the azimuthal asymmetry has been pa- 
rameterized following a simple hydrodynamically motivated blast wave model, 
described in [19], [43], and generalized in [42] to also describe v\ : 



2?r 

/ rf0 6 cos(n0 fe )/ n (a)Ki(/3)[l + 2s n cos(n0 6 )] 

Vn(PT) = (8) 

/ # 6 io(a)-^i(/3)[l + 2s„ cos(n0 fe )] 



o 



where the harmonic n can be either 1 or 2, I , I n , and K\ are the modi- 
fied Bessel functions, 0^ = — ^rp, = (pr/Tf) sinh[p(0 b )], f3{4>b) = 
(mr/Tf) cosh[p(0b)], s n is the surface emission parameter, and Tf is the freeze- 
out temperature. The azimuthal flow rapidity is given as p(<f>b) = Po+Pa cos(n0b) 
with p is the mean transverse expansion rapidity (vo = tanh[p ]) and p a is the 
amplitude of its azimuthal variation, respectively. Further details are given in 
[42]. 

Without regard to the physical interpretation of the parameters of this model 
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it can be used to provide a convenient parameterization of the dependence 
of the flow. For this purpose Eq. 8 can be further simplified by setting p a = 0. 
The integrals can be solved to reduce the expression to two parameters, a n 
and b n . The expressions for the p T dependence of the directed and elliptic flow 
are then: 



/ x h{biP T ) (a s 



The measured photon v n (pr) are described well by these expressions (see 
Fig. 4c, d). The extracted a n and b n were constrained to have a smooth cen- 
trality dependence. 

The simulated n° pr distributions have been parameterized as: 



d 2 N 1 d*N ( p' \ n 

-E-t^Pt = Cp T \ -7- (11; 



N Event dydpT N Event dp 3 yp'o+Pr, 

with parameters C',p' ,n' taken from the measured ir° results ([26]). Jetset 
7.4 was used to generate the 7r°-decay photons which were then filtered in the 
simulations with the detector acceptance and efficiency 

In the simulations, the 7r°'s were generated with an azimuthal asymmetry with 
Pt dependent v n parameterized by Eqs. 9 and 10. As an initial ansatz, the 
measured photon v n = g(pr) for each centrality was used for the v n = f(pr) of 
the simulated pions. The output v n = g'{pr) of the simulated decay photons 
was compared to the measured photon dependence and used to adjust the 7r° 
v n for the next iteration. The procedure was iterated until the simulated and 
measured photon results were in agreement, which typically required two to 
three iterations. The a n and b n coefficients obtained from the fit to the photon 
results and the coefficients extracted for the 7r°'s by the iteration procedure 
are summarized in Table 2. The ratio of the flow coefficients extracted for 
the simulated decay photons to the input pion flow coefficients provided the 
correction factors k n {pr) = v Z(.Pt) / VniPr) used to extract the 7r° v n from 
the measured photon v n for p T > GeV/c. The systematic errors on the k n 
were determined from the uncertainties A of the fit parameters: a n ± Aa n , 
b n ± Ab n , p' ± Ap' and n' ± An'. Within errors, the final results were found 
to be consistent with an analysis in which a simple linear pt dependence of 
the v^(pt) was assumed. 
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Table 2 

The blast wave parameters a n and b n of Eqs. 9 and 10 for various centralities 
obtained by fits to the WA98 inclusive photon data (left) and the corresponding 
parameters extracted for neutral pions (right) from simulations that reproduce the 
inclusive photon results. 



Fig. 5a,b shows k 2 , the ratio of the flow coefficients of photons and pions, 
obtained from the simulations. The k\ ratio is found to be independent of 
centrality and pr with a value of (ki) = 1.075 ± 0.003 with a systematic error 
less than ±12%. Since v\ is small, k\ « 1. The k 2 ratios are also found to be 
independent of centrality with (/c 2 ) = 1.215 ±0.002 and a systematic error less 
than ±8%, but strongly dependent on p?. The observation that the k n ratios 
are greater than unity can be understood as a simple effect of the n° decay 
and the fact that the v n increase with p T . Comparing photons and 7r°'s at a 
given p T , the decay photons will have been produced from 7r°'s with larger p T 
and hence larger v n , giving k n > 1 (see Fig. 5a and b). The k n results shown 
in Fig. 5a and b have been used to estimate the neutral pion flow values from 
the measured photon flow. 

The 7T° t> 2 extracted from the measured WA98 photon v 2 as a function of 
centrality are compared to the results from NA49 for charged pions [42] in 
Fig. 6. For this comparison the WA98 results for photons with > 0.18 
GeV/c were corrected to the expectation for the pions without pr threshold 
(Pt > GeV/c) as for NA49 results with the correction factor k® = 1.30±0.036 
and k\ = 1.59 ± 0.003. The systematic errors on the WA98 points including 
the errors of the measured photon v n values and the additional error of the 
correction are less than ±20%, except for the lowest pt points where the 
upper systematic errors increase (Fig. 8). The systematic errors on the NA49 
data points vary from 13% for the most peripheral to 80% for the most central 
bin. The two measurements are seen to be in agreement. 
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The rapidity dependence of the pion flow coefficients is shown in Fig. 7 for 
two centrality selections. The WA98 photon v n with > 0.18 GeV/c were 
corrected to the expectation for 7r°'s without pt threshold using the correction 
factors k® given above. The NA49 results for pions with p T > GeV/c with 
appropriate corrections ([42]) are also shown and seen to be in agreement. 
Similar agreement between the two experiments is seen in the transverse mo- 
mentum dependence of the v 2 shown in Fig. 8. 

Fig. 9 shows the WA98 ir° v 2 coefficient deduced from the measured photon 
v 2 (corrected by k®) together with a compilation of results [42] from other 
experiments. Results are shown from E877 [8,12], CERES [44], NA49 [42], 
PHOBOS [45], PHENIX [46], and STAR [47], The WA98 result follows the 
general trend of the smooth increase of the elliptic flow with increasing beam 
energy. It should be kept in mind that the pt and rapidity coverages as well 
as centrality selections differ for the various experiments. 

Since the strength of the elliptic flow should be proportional to the initial 
eccentricity of the collision zone, it is useful to normalize the measured v 2 to 
the eccentricity of the reaction geometry when investigating the systematics 
of the elliptical flow [48]. The initial spatial eccentricity was calculated within 
a Glauber model calculation [37] as: 

_ < y 2 > - < x 2 > 
< y 2 > + < x 2 > 

where x and y are the participant nucleon coordinates in the plane perpendic- 
ular to the beam and x denotes the in-plane direction. 

It is of interest to plot the quantity v 2 /e versus the particle density as estimated 
by the dN/dy of charged particles divided by the area of the overlap region 
S [8,42]. Neglecting the weak incident energy and centrality dependence of 
the average pt, the scaled particle density is proportional to the initial energy 
density in the Bjorken estimate [51]. Fig. 10 shows the WA98 result for v 2 /e 
for neutral pions together with results for charged pions from NA49 [42] and 
results for charged particles from STAR [49,50] and E877 [8]. The RHIC data 
have been corrected for their p T cutoff, the errors are statistical. For the WA98 
points the dN/dy are taken from WA98 measurements [37]. 

The WA98 results confirm the previous observation [8,42] of a universal depen- 
dence of v 2 /e on the particle density. It should be noted that results at different 
incident energies but similar average particle density correspond to dramati- 
cally different collision geometries. Generally, a central collision at low energy 
will have a particle density similar to a more peripheral collision at higher 
energy. Thus, the universal dependence demonstrates that the scaled v 2 de- 
pends only on the initial energy density, or initial pressure, rather than being 
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dependent on the particle density times pathlength, which would be the ex- 
pectation if only partial thermalization was attained. The observed increase of 
the scaled elliptic flow with increasing particle density is qualitatively similar 
to that seen in a recent hydrodynamical model study [52] in which the v 2 value 
for fixed impact parameter increases smoothly with increasing dN c h/dy corre- 
sponding to the increase in pressure. That systematic study of hadron spectra 
and flow results indicated that existing data are best reproduced by hydrody- 
namic model calculations with an equation of state which include a transition 
to a Quark Gluon Plasma Phase at a critical temperature of T c = 165 MeV 
with a latent heat of 800 MeV [52]. 



7 Summary 

Directed and elliptic flow of photons have been measured in 158 A GeV Pb+Pb 
collisions using the LEDA electromagnetic calorimeter and the Plastic Ball 
detector of the WA98 experiment. The conventional reaction plane method 
for directed and elliptic flow and the pair-correlation method for elliptic flow 
have been used in the analysis. The elliptic flow values obtained by the two 
methods are consistent. 

The centrality and pt dependences of the photon directed and elliptic flow 
coefficients were presented. In-plane elliptic flow is observed. Both flow values 
increase for more peripheral collisions and with increasing px- 

For comparison with charged pion measurements, the neutral pion flow coef- 
ficients have been extracted from the measured photon flow coefficients using 
Monte Carlo simulations. The Monte Carlo simulations demonstrate that the 
average photon flow coefficients are ps7-21% greater than the parent 7r° flow 
coefficients. The extracted neutral pion flow results are compatible with the 
NA49 charged pion flow and with the general trend of the elliptic flow be- 
haviour as a function of beam energy. The universal dependence of the elliptic 
flow, scaled by the eccentricity of the initial nuclear overlap, on particle density 
is confirmed. 
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Fig. 1. Left: The measured photon azimuthal correlation functions for different 
centralities. The solid lines show fits to equation 3. Right: Subevent correlation 
functions of particles measured in the Plastic Ball for different centralities. The 
filled circles show the measured correlation functions. The open circles show the 
results for mixed subevents. 
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Fig. 2. (a) The resolution correction factors obtained from subevent correlation func- 
tions of particles measured in the target fragmentation region as a function of the 
number of participants. The filled circles show the values for directed flow, the open 
circles for elliptic flow, (b) The elliptic flow V2 of photons extracted by the reaction 
plane method (solid circles, pt > 0.18 GeV/c) and by the correlation method (open 
circles, p T i > 0.18 GeV/c, p T 2 > 0.18 GeV/c), integrated over p T > 0.18 GeV/c 
and y =2.3-2.9 as a function of the number of participants. 
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Fig. 3. The elliptic flow V2 of photons extracted by the reaction plane method (solid 
circles, px > 0.18 GeV/c) and by the correlation method (open circles, pri > 0.18 
GeV/c, pT2 > 0.18 GeV/c), integrated over j/ =2.3-2.9 for three centrality selec- 
tions. 
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Fig. 4. Directed flow (vi) (a) and elliptic flow (t^) (b) of photons integrated over 
Pt > 0.18 GeV/c and y = 2.3 — 2.9 as a function of the number of participants. 
Directed flow (vi) (c) and elliptic flow (t^) (d) of photons integrated over y =2.3-2.9 
as a function of pt for various centralities. Solid circles show results for semi-central 
events (47-83%), open circles show results for semi-central events (13-47%), solid 
squares show results for central events (0-13%). Solid lines are fitted results to a 
blast wave model functional form. 
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Fig. 5. The ratio of photon elliptic flow to parent pion elliptic flow coefficients, (a) 
&2 = vj/v^ for min. bias, (b) A?2 = vj/v^ for different centrality bins. 
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Fig. 6. The centrality dependence of pion elliptic flow V2 extracted from the mea- 
sured photon flow from WA98 (solid circles, —0.6 < ycu < 0), compared to pion 
flow results of NA49 (open circles, < ycu < 2.1, [42]). 
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Fig. 7. The rapidity dependence of pion flow extracted from the measured photon 
flow from WA98 (circles) compared to pion flow results of NA49 (stars, [42]). The 
results are shown for central on the left panel (WA98: 0-13%, NA49: 0-12.5%) and 
semi-central on the right panel (WA98: 13-47%, NA49: 12.5-33.5%) selections for 
the directed flow (a,b) and elliptic flow (c,d). Systematic error bands are shown on 
(c) for NA49 (dashed curves) and for WA98 (solid curves). Open points are reflected 
at mid-rapidity. The dashed lines indicating the rapidity dependence are only meant 
to guide the eye. 
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Fig. 8. The transverse momentum dependence of the tt elliptic flow extracted from 
the measured photon flow from WA98 (solid circles) compared to ir + elliptic flow 
results of NA49 (open circles, [42]). The results are shown for semi-central selections 
(top: WA98: 13-47%, -0.6 < y C u < 0, NA49: 12.5-33.5%, < y C u < 0.8) and 
for central selections (Bottom: WA98: 0-13%, -0.6 < y C M < 0, NA49: 0-12.5%, 
< Vcm < 0.8). Systematic error bands are shown for NA49 (dashed curves) and 
for WA98 (solid curves). 
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Fig. 9. Comparison of the WA98 p-p-integrated elliptic flow results at \fs^N = 17 
GeV (13-47%, —0.6 < ycM < 0) with results from other experiments for differ- 
ent collision energies. The WA98 result is the 7T° elliptic flow extracted from the 
measured photon elliptic flow, as described in the text. The total statistical and 
systematic error is shown on the WA98 point. 



27 



>0.25 



0.2 



0.15 



0.1 



0.05 











• WA98, E |ab /A=158 GeV, y CM : -0.6 - 0.0 
O NA49, E |ab /A=158 GeV, y CM : 0.0 - 0.8 




□ E877, E |ab /A=11.8GeV 
A NA49, E |ab /A=40 GeV 
STAR, VS NN =130GeV 
* STAR Prelim., VS NN =200 GeV 



10 



15 



20 



25 



30 



35 



Pch =(1/S)dN ch /dy(fm- 2 ) 



Fig. 10. The scaled elliptic flow i^/e as a function of particle density. Results are 
shown for the WA98 pr-integrated V2 values at ^/snn = 17 GeV ( —0.6 < ycM < 0). 
The WA98 results are compared to results from the E877, NA49, and STAR exper- 
iments. The WA98 result is the 7r° elliptic flow extracted from the measured photon 
elliptic flow, as described in the text. 
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